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Thermoelectric properties of the system La2Ni04+<5 have been studied ab initio. Large Seebeck 
coefficient values are predicted for the parent compound, and to some extent remain in the hole- 
doped metallic phase, accompanied of an increase in the conductivity. This system, due to its 
layered structure would be a suitable candidate for an improvement of its thermoelectric figure of 
merit by nanostructurization in thin films, that has already been shown to increase the electrical 
conductivity (a). Our calculations show that in the region around La2NiO4.05 the system has a large 
thermopower at high temperatures and also a substantially increased a. Films grown with this low- 
doping concentration will show an optimal relationship between thermopower and a. This result 
is obtained for various exchange-correlation schemes (correlated, uncorrelated and parameter-free) 
that we use to analyze the electronic structure of the hole-doped compound. 



I. INTRODUCTION 

In recent times, the studies on the thermoelectric 
properties of various compounds both using exper- 
iments, theoretical and computational tools has fo- 
cused in the search for new functional materials that 
could be used in the fabrication of energy-recovery 
devices based on a thermoelectric effect. Oxides are 
promising in principle because of their chemical sta- 
bility and possibility of usage in high temperature 
devices. Particularly, misfit layered cobalt ates have 
drawn significant attention^ - — Large thermopower 
together with metallic conductivities place that fam- 
ily as one of the most relevant and thoroughly stud- 
ied within the oxides. One electronic structure rea- 
son behind their interesting thermoelectric proper- 
ties is the existence of two electronic systems within 
those cobaltates. 6 Due to the trigonally distorted oc- 
tahedral environment, the \>2g bands are split be- 
tween a highly mobile wide band that drives metal- 
lic conduction (of e* character) and a more localized 
band (with ai p parentage), narrow in nature that 
provides a large Seebeck coefficient i 

The search for new oxide compounds that could 
present a similar thermoelectric response is an inter- 
esting issue to be pursued. In this paper we will ar- 
gue that a possible candidate family for that sake are 
a particular kind of layered nickelates: hole-doped 
La2Ni04. This compound and others related to it 
have drawn the attention of the community for var- 
ious reasons: i) La2Ni04 is isostructural with the 



parent compound of the high-T c cuprates: 8 MO2 
planes with a square planar coordination around the 
metal cations form the building blocks of this lay- 
ered material, with La acting as spacer between the 
layers, ii) an interesting temperature-induced metal- 
insulator transition (MIT) occurs in it^ iii) it is pos- 
sible to lower further the valence of the Ni cation 
in a layered structure to compare it with supercon- 
ducting cuprates, not only in terms of structure but 
also with the same electron count iSr— Back in the 
eighties, the MIT was explained by Goodenough and 
Ramasesha 18 as occurring in a system where two dif- 
ferent d-states coexist around the Fermi level, hav- 
ing clearly distinct band widths. This is reminiscent 
of the above mentioned cobaltates, but in this case 
two different e^ bands exist, a more localized d z 2 
and a more itinerant d x i_ y i band. These yield rea- 
sonable conductivities, and also due to the existence 
of a more localized set of narrow bands, interesting 
values of the thermopower could emerge. 

Experimental studies of the Seebeck coefficient 
in La2Ni04 show promising values exceeding 100 
/iV/K at room temperature^ (the typical minimum 
value for a material to be considered a high per- 
formance thermoelectric) even though most reports 
with oxygen-rich samples show smaller values (we 
will argue below what the reason for this is). The 
electrical conductivity is larger in the MO2 planes, 20 
which is the expected behavior for such a layered 
structure. To obtain an estimate of the thermoelec- 
tric performance of a material, it is common to use 



the so-called figure of merit, which is defined as the 
dimensionless quantity zT=crTS 2 / n. zT > 1 is re- 
quired for applications. 21 With the best values ob- 
tained experimentally for S ~100 /iV/KjiS p ~ 5 
mf].cm2^ and k ~ 7.5 W/m.K, 23 a zT at room 
temperature of about 0.01 can be estimated, which 
is comparable to other promising compounds such as 
CrN j 24 i 25 particularly so because of the large room 
for improvement, some of which we will discuss in 
this paper. In particular, the layered structure of 
the compound would be prone to yield better con- 
ductivity when grown in the form of thin films, and 
this has been shown in the past* 2 ^ A thin film geom- 
etry would hamper thermal conductivity along the 
c-axis, which would eventually improve the thermo- 
electric figure of merit as long as the thermopower 
and the mobility are not reduced.— Also, oxygen 
excess can lead to hole doping and a corresponding 
increase in electrical conductivity. A vast amount of 
experimental information exists on the thermoelec- 
tric and transport properties of La2Ni04. However, 
an exact account of the oxygen content is necessary 
to obtain good systematics, and this is not easily 
found in literature, probably due to complications 
in the correct determination and stabilization of a 
particular oxygen content. 

In this paper we will explore a possible path to- 
wards an improved thermoelectric response via oxy- 
gen doping. We will discuss in particular how mov- 
ing the system into the more itinerant limit by 
hole doping can reduce further the importance of 
the phononic part of the thermal conductivity (that 
works against a good thermoelectric performance) 
relative to the electronic part, and also improving 
the figure of merit if a decrease in the thermopower 
is prevented. Concerning the magnetic order, it has 
been shown that La2Ni04 is an antiferromagnet with 
an in-plane ordering such that an antiferromagnetic 
(AF) interaction between nearest neighbor Ni atoms 
is stabilized. The commensurability of the order- 
ing has been put into question, 27 but careful studies 
with respect to the oxygen content show how the 
Neel temperature is significant even at values of S ~ 
0.1^ and could vanish at S ~ 0.14J^ We will focus 
our calculations in that doping level (5 < 0.15), close 
to the localized electron limit where magnetic order 
is present at low temperature and our DFT calcu- 
lations, that assume a magnetically ordered phase 
at zero temperature can describe the system more 
accurately. All calculations presented here will have 



that type of ordering imposed. 

The paper is organized as follows: Section lU will 
present the calculation details, Section mil will de- 
scribe the electronic structure, analyze the thermo- 
electric properties calculated for various oxygen con- 
tents, and finally we will summarize the main con- 
clusions of the work in Section HVl 



II. COMPUTATIONAL PROCEDURES 

Our electronic structure calculations were per- 
formed within density functional theor y 29 ! 30 using 
the all-electron, full potential code wien2k 31 based 
on the augmented plane wave plus local orbitals 
(APW+lo) basis set. 32 We will present calculations 
of the electronic structure and transport properties 
obtained with various exchange-correlation function- 
al and for several doping levels. For the exchange- 
correlation functional in the calculations at U= 0, 
we have used the Perdew-Burke-Ernzerhof version 
of the generalized gradient approximation 33 (GGA). 

To deal with strong correlation effects that are 
widely acknowledged to play an important role 
in this type of nickelates, we apply the LDA+/7 
schem o 34 i 35 that incorporates an on-site Coulomb 
repulsion U and Hund's rule coupling strength J# 
for the Ni 3d states. For the uncorrelated part of 
the exchange-correlation functional we used the local 
density approximation (LDA) 36 in that case. The 
LDA+£7 scheme utilized is the so-called "fully lo- 
calized limit" 37 (FLL). A description of the results 
obtained at different values of U (in a reasonably 
broad range 4.5-8.5 eV for the Ni cations) was per- 
formed, but results presented here will be for U= 7 
eV. No significant changes in the main conclusions 
of the paper are found within that range of U values. 
For a discussion of the effects of a different choice of 
U and LDA+U method in similar layered nickelates, 
see Ref. [38|. The value chosen for the on-site Hund's 
rule strength is J = 0.68 eV and is kept fixed. 

Also, the recently developed Tran-Blaha modified 
Becke- Johnson (TB-mB JLDA) potential was used^ 
This potential has been shown to provide an ac- 
curate account of the electronic structure of corre- 
lated compounds 40 using a parameter-free descrip- 
tion, also yielding accurate band gaps for most semi- 
conductors. 

All calculations were fully converged with respect 
to all the parameters used. In particular, we used 
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K rnt K rnax = 7.0, a k-mesh of 10 x 10 x 4, and muffin- 
tin radii of 2.35 a.u. for La, 1.97 a.u. for Ni and 1.75 
a.u. for O. 

The transport properties were calculated using 
a semiclassical solution based on the Boltzmann 
transport theory within the constant scattering time 
approximation by means of the BoltzTrap code, 
that uses the energy eigenvalues calculated by the 
wien2k code. 41 In this case, denser k-meshes are 
required, in our case up to 40 x 40 x 15 to reach 
convergence. The constant scattering time approxi- 
mation assumes the relaxation time r(e) as energy- 
independent. This results in an expression of both 
the thermopower and the thermoelectric figure of 
merit with no dependence on r (it can be directly 
obtained from the band structure without any as- 
sumed parameters). This approximation has been 
used succesfully to describe several thermoelectric 
materials 42r— 

To simulate small doping values we have used the 
virtual crystal approximation (VCA)4£ This con- 
sists in modifying the total number of electrons in 
the system (also the total atomic number for charge 
neutrality) to simulate a particular doping level. 



III. RESULTS 
A. Basic description of the electronic structure 

A simple ionic picture of La2Ni04 gives a formal 
valence Ni 2+ :d 8 , which is always a S= 1 ion. Typi- 
cally, for this type of low- valence nickelate, with the 
Ni 2+ cations in a square planar environment, one 
would expect the xz/yz doublet fully occupied, and 
a higher position of the xy singlet competing in oc- 
cupation with the e g bands (which are also split and 
will have significantly different band widths). The 
typical situation would be a lower xy band fully oc- 
cupied and two e^ unoccupied bands, being the band 
gap 4 ^ of a distinct d-d character. 

Figure [1] shows the calculated density of states 
(DOS) for La 2 Ni0 4 in an AF solution, both for GGA 
and a TB-mBJLDA scheme together with the GGA 
band structure in the lower panel. Similar DOS 
curves are obtained with the various functionals, the 
size of the gap being the main difference, and also 
some differences exist in the particular bandwidths 
and position of the bands, but even calculations with 
U= yield an insulating solution as long as in-plane 
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FIG. 1: (Color online.) DOS of the Ni e g levels in 
La 2 Ni0 4 ordered AF within GGA and TB-mBJLDA. 
Observe the more localized d z 2 band and the broader 
d x 2- y 2, particularly noticeable in the GGA calculation. 
Lower panel shows the band structure in the GGA cal- 
culation. The t2 9 bands are in the same energy region 
as the majority d z 2 band. 



AF order is set. In the plots, we only show the e^ 
states of the Ni cations since these have the domi- 
nant contribution above and immediately below the 
Fermi level (and we want to see them with some 
degree of clarity). Furthermore, due to the square 
planar oxygen environment, the e g states are sig- 
nificantly split, and differ substantially not only in 
position but in width. The d x i_ y i band is higher 
in energy due to crystal-field splitting caused by the 
absence of apical oxygens around the Ni cations, but 
it is slightly wider, whereas the d z i band is lower but 
it is composed of more localized electrons (it is nar- 
rower). One can see this very clearly in the band 
structure figure. The d z i bands appear both in the 
GGA and TB-mBJLDA calculations much narrower, 
with less than 1 eV band width, being the corre- 
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sponding width of the d x 2_ y 2 band at least 2 eV. In 
the band structure figure (lower panel), we can also 
see that not only the e g levels are in the vicinity 
of the Fermi level but also the t2 g states are placed 
somewhat close to it, and are actually in the same 
enery region as the majority-spin d z 2 band. 

This is a narrow-gap oxide (of only 70 meV 48 ac- 
cording to resistivity measurements), but calcula- 
tions show the gap is substantial (0.6 eV for GGA, 
1 eV for LDA+Z7 at U= 7 eV and 2 eV for TB- 
mBJLDA). Other narrow-gap oxides show this type 
of band gap overestimation coming out from cal- 
culations (CrN, 40 SreCosOi^) and still the calcu- 
lations of the transport properties agree well with 
experiments. One possibility for this could be the 
fact that the low gap value obtained from resistiv- 
ity measurements comes from shallow states in the 
gap, being the actual band gap much larger, as hap- 
pens in CrN. 40 In the case of La2Ni04, a small band 
gap is consistent with the temperature-induced MIT 
observed. The MIT in this system has been de- 
scribed as a temperature-induced broadening of the 
d x 2_ y 2 band that eventually overlaps the occupied 
with the unoccupied band of opposite spin yielding 
metallic behavior at high temperature in the parent 
compound. 9 Our calculations confirm the descrip- 
tion of the electronic structure that Goodenough and 
Ramasesha have done in the past, 18 although we see 
that the t2 g bands are located together with the ma- 
jority d z 2 band in the vicinity of the Fermi level. The 
GGA DOS we show would anticipate that a sufficient 
temperature-induced broadening of the x 2 — y 2 band 
would place the unoccupied band edge below that of 
the unoccupied z 2 band, which is significantly nar- 
row. Something similar can be expected for the TB- 
mBJLDA calculation, but the largely overestimated 
gap makes the analysis a bit more adventurous. 



B. Tuning the electronic structure via 
hole-doping 

As mentioned before, if somehow one can com- 
plement metallic conductivity with the maintenance 
of a large thermopower and a reduced thermal con- 
ductivity, the material could then be comparable to 
other oxides in thermoelectric performance^ Let us 
see how this can be done in this particular case. As 
can be seen, the DOS rises rapidly right below the 
top of the valence band (d x 2_ y 2 -character) which is 



generally favorable for the thermopower. A similar 
heavy-electron band appears somewhat also in the 
electron-doping region just above the bottom of the 
conduction band (in that case, of d z 2 parentage) al- 
though with less dispersive character. Hence, the 
more localized d z 2 band should contribute substan- 
tially to the Seebeck coefficient in the hole-doping 
region as well, but its DOS peak is very far from the 
Fermi level to be attained by a relatively small hole- 
doping level. This would instead move the chemi- 
cal potential into the occupied x 2 — y 2 band with a 
noticeable increase in electrical conductivity. This 
band could also be sufficiently narrow to retain a 
substantial value of the thermopower if the chemi- 
cal potential lies in a region inside the valence band 
comprising a large derivative of the density of states. 

In order to see this in detail, we can start by 
analyzing the electronic transport properties of the 
parent compound. We have calculated these for 
La2Ni04 using various exchange-correlation poten- 
tials. Our calculations were carried out using the 
semiclassical approach, within the constant scatter- 
ing time approximation. Yet, the value of the ther- 
moelectric figure of merit zT= T<rS 2 /ft is indepen- 
dent of the scattering time chosen and so are the 
thermopower and ct/k, as long as the constant scat- 
tering time approximation is retained. In Fig. [2] we 
present the thermoelectric figure of merit as a func- 
tion of the chemical potential calculated at 400 K 
using different exchange-correlation functionals. We 
study the electronic part only, we do not consider 
phonon terms in the thermal conductivity, which 
might be substantial in particular for the localized 
electron limit. Thus, this electronic-only zT can be 
considered as a theoretical upper limit of the figure 
of merit for this system. We will see below how we 
can make an estimate of these effects based on the 
experimental data for a and n from Refs. Hq|, [H and 

Ef 

Figure [2] shows a reasonably large value of the 
thermoelectric figure of merit for the stoichiomet- 
ric compound (chemical potential at the top of the 
valence band), which is given by the large positive 
Seebeck coefficient, in qualitative agreement (see the 
discussion below) with experimental observations at 
low-J (oxygen excess) valuesJ^ However, this is ex- 
pected to be hampered by the large thermal conduc- 
tivity due to phonons (ftph), which would substan- 
tially surpass its electronic counterpart (ft e j), which 
is the only component we are analyzing with these 
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FIG. 2: Thermoelectric figure of merit (electronic part only) of La2Ni04 AF within GGA, the LDA+?7 approach and 
TB-mBJLDA as a function of chemical potential represents the Fermi level). Observe a peak in the figure of 
merit for small hole-doping levels. zT becomes zero at about -0.3 eV when the Seebeck coefficient changes sign from 
positive to negative, a similar feature happens in the electron-doped region, S goes from negative to positive when 
the system is further electron doped. 
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FIG. 3: (Color online.) Calculated thermopower at T= 
400 K as a function of hole doping in La2Ni04+<5 AF 
within GGA, the LDA+[7 approach and TB-mBJLDA . 



calculations. In any case, moving towards the itiner- 
ant electron limit, when K e \ starts being comparable 
to K p h, these results become more significant. We 
observe an interesting peak in the figure of merit in 
the hole-doping region (oxygen excess, positive val- 
ues of 5, negative chemical potential). This peak 
is present independently of the exchange-correlation 
functional used, but the characteristics (height and 
position in chemical potential) vary depending on 
the method used. 



In this paper we try to describe a physical phe- 
nomenon that occurs in a range of physical parame- 
ters (in this case the oxygen excess) in La2Ni04 + 5 
but the quantitative values of which are very dif- 
ficult to obtain ab initio. The peak is related to 
the x 2 — y 2 bands, as we can see comparing Fig. [2] 
with the DOS plots in Fig. [1] (the energy window 
of the large peak in zT lies inside that band). We 
also observe that other peaks in zT appear at higher 
hole-doping concentrations (more negative chemical 
potential), related in that case with the d z i and 
bands. Even though the d z i band is certainly nar- 
rower, it does not provide a sufficiently high zT com- 
pared to x 2 — y 2 (a multipeak structure is predicted 
by the TB-mBJLDA functional) because it is merged 
in the same energy region as the t2 g bands. If some- 
how band engineering could be done to bring the 
d z 2 band closer to the Fermi level, the possible ther- 
mopower obtained could be enhanced taking benefit 
of the localized nature of the band, and at a signif- 
icantly smaller doping level. The realization of this 
possibility looks difficult. It can also be oberved 
how zT becomes zero at about -0.3 eV simultane- 
ously with the Seebeck coefficient sign change from 
positive at low doping to negative at higher doping 
levels. A similar feature happens in the electron- 
doped region, S goes from negative to positive when 
the system is further electron doped (about 0.2 eV 
above the bottom of the conduction band). 

In order to characterize in more detail that first 
peak in zT at 5 ^ 0, and see how the amount of 
oxygen excess (hole-doping) relates to a particu- 
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lar change in chemical potential in the system we 
present in Fig. [3] the calculated thermopower of 
the hole-doped system at 400 K for all the func- 
tional used. Since the behavior of a j n is relatively 
monotonous both with temperature and doping, the 
behavior of the thermopower (its peaks and sign 
changes) gives the broad picture of how the figure of 
merit itself evolves. The dependence of the Seebeck 
coefficient with both the doping level and the S value 
is shown. A peak can be observed in all the curves 
although it is more pronounced for LDA+Z7 and be- 
comes broadened within TBmBJLDA. The height 
of the peak is related to the width of the x 2 — y 2 
band, which is narrowed by the LDA+/7 method, 
and broadened by the TB-mBJLDA scheme. An in- 
termediate situation is given by uncorrelated GGA. 
The doping level necessary to reach it depends on 
the relative position of that band with respect to 
the Fermi level in the calculations, GGA placing it 
further away than the other methods (see Fig. [TJ 
top panel). It can be seen that for LDA+/7 the 
peak is located at 5~ 0.03, whereas for GGA and 
TB-mBJLDA it appears at higher doping levels ~ 
0.06 and 0.07, respectively. In addition, the above 
described sign change of the Seebeck coefficient in 
the hole-doped region can be seen for both GGA and 
LDA+/7 at S close to 0.15 (larger, beyond the values 
of doping we are studying here, for TB-mBJLDA). 
This has been observed experimentally 19 for values 
around 5 ~ 0.20, in good agreement with the result 
obtained from our calculations. All these calcula- 
tions were obtained for the parent undoped com- 
pound La2Ni04 as a function of the chemical poten- 
tial, what could be called a rigid-band approach. We 
will argue below how the picture is modified when 
VCA is applied. 



C. VCA calculations 

Once we have localized the position of the peak 
as a function of doping, we will study the evolution 
with hole-doping of the Seebeck coefficient (com- 
paring it with available experimental data) and the 
other transport properties that can be calculated. 
In order to do that, we have performed VCA calcu- 
lations (modifying the total number of electrons in 
the system to simulate the dopant concentration). 
Since we will be dealing with small doping levels (5 
< 0.15), we will consider the ground state AF con- 
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FIG. 4: (Color online). Thermopower of La 2 Ni0 4 AF 
within GGA, the LDA+J7 approach and TB-mBJLDA 
and experimental values as a function of temperature for 
various values of the hole doping concentration S simu- 
lated using the virtual crystal approximation. Observe 
the different tendency with temperature for optimum- 
doping values (increasing at high temperature) and 
higher values of doping (oscillating with temperature and 
flat at the high temperature end). These tendencies are 
consistent with experiments (shown in the lower right 
panel) . 



figuration and also the structural properties of the 
parent compound. With this, we neglect the effects 
of the known orthorhombic distortio n 50 ^ 1 that oc- 
curs in this nickelate and would be a second order 
effect for the properties we try to estimate. The 
shape of the curves in Fig. [2] will not remain the 
same when VCA is applied because the effect of dop- 
ing is not just a rigid band shift, but the existence 
of an optimum zT with doping remains true. Also, 
the more pronounced maximum in figure of merit at 
negative values of the chemical potential occurs for 
the LDA+Z7 calculations, as we have seen in Figs. [2] 
and [3] and will further confirm below. 

Figure |4] shows the thermopower at different dop- 
ing levels calculated within the various exchange- 
correlation potentials we are using, and an addi- 
tional figure with some experimental data from the 
literature. We can observe large values of the See- 
beck coefficient for the most insulating samples, 
where the thermopower is above 200 /iV/K for the 
whole temperature range. The large Seebeck coeffi- 
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FIG. 5: (Color online). Thermoelectric figure of merit (electronic only) of La2Ni04 AF within GGA, the LDA+U 
approach and TB-mBJLDA as a function of temperature for various values of the hole doping concentration S 
calculated using the virtual crystal approximation. Observe the large high temperature figure of merit attainable at 
small doping values. 



cient at low doping, 19 is well predicted by all func- 
tionals. In every case, the maximum in zT occurs in 
the doping region of our interest (S < 0.15) for val- 
ues that coincide roughly with the maxima in ther- 
mopower vs. 5 at 400 K we analyzed in Fig. [3] (below 
5 ~ 0.05 within LDA+U and at slightly higher values 
for GGA and TB-mBJLDA). For these "optimum" 
doping values the thermopower increases at high 
temperatures while for the other 5 values it shows 
an oscillating behavior and a flattening at higher 
temperatures, results that are compatible with most 
experiments available. The oscillating behavior of 
the thermopower with temperature that has been 
found in literature is particularly well described by a 
GGA calculation at x ~ 0.05 and also TB-mBJLDA 
at 5 ~ 0.10. According to our calculations, exper- 
iments showing thermopower oscillating with tem- 
perature (showing a maximum at intermediate tem- 
peratures) have been obtained analyzing samples in 
that doping region (hole-doped, non-stoichiometric, 
even if some of them do not explicitly mention the 
oxygen content) j^ 2 - - — However, an accurate place- 
ment of the maximum in the Seebeck coefficient as a 
function of temperature is not given by our calcula- 
tions, that type of agreement being similar to what 
has been found in other correlated oxides. 49 Zhou 
et al. 52 have shown a slight increase in S at high 
temperature experimentally, but the exact 5 value 
of their samples is not provided; from our calcula- 
tions we can argue that those samples would have 
a substantial S- value on the order of 5 ~ 0.15-0.2. 
The best quantitative agreement between our calcu- 
lations and those reported from experiments at low- 



doping levels can be obtained with the TB-mBJLDA 
functional. This was also found for the thermoelec- 
tric properties of CrN, when compared with exper- 
iments and with LDA+/7 calculations. 40 This func- 
tional provides a purely ab initio parameter-free de- 
scription of the system, not depending on the choice 
of an arbitrary parameter. 

In Fig. [5] we present the evolution with tempera- 
ture of the thermoelectric figure of merit (electronic- 
only) for various doping levels. The figure com- 
plements the description given by the thermopower 
analysis alone, with similar conclusions. LDA+U 
values for 5= 0.05 would be very promising in terms 
of the figure of merit obtained, even surpassing that 
of the undoped compound at high T. Also the values 
predicted by the parameter-free TB-mBJLDA func- 
tional would exceed zT= 0.5 at high temperatures 
for 5= 0.05. Uncorrelated GGA predicts an opti- 
mum hole-doping content at 5= 0.10 with large fig- 
ure of merit at high temperature. LDA+/7 also gives 
a very good response for 5= 0.10 and 0.15, especially 
due to its high temperature dependence. From Figs. 
[2] and [3j focused at 400 K and carried out with a 
rigid-band approach, we expected a sharper behav- 
ior in S- values, but the VCA results show a broader 
range of efficient thermoelectric response being pre- 
dicted. We have already discussed that the large zT 
of the 5= case would need to be reduced by intro- 
ducing a substantial n p h- This would be particularly 
important at 5= 0. Being the system an insulator, 
K p h should be much higher than its electronic coun- 
terpart. Our calculations at S ^ should be more 
realistic, being the system approaching the metallic 
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FIG. 6: Thermoelectric figure of merit of La2Ni04 AF as a function of temperature for various values of the hole 
doping concentration 5 within GGA, the LDA+t/ approach and TB-mBJLDA estimated using the experimental 
values of k and a 



limit, specially if the thermal conductivity due to 
phonons can be further reduced by nanostructuring. 
This has been shown in the past to be a key ingredi- 
ent in the design of high performance thermoelectric 
devices. 55 In that case, i^d could become more im- 
portant compared to K p h, as discussed above. Let us 
give some rough numbers. According to our calcu- 
lations (within GGA, not shown), the conductivity 
increases by a factor 20 at room temperature for 5= 
0.05 with respect to the stoichiometric compound, 
and by a factor 80 for 5= 0.15, explainable by the 
increase in carrier concentration that is also given in 
Fig. [21 Also, if one can make thin films of the appro- 
priate oxygen excess concentration, the conductivity 
can be made even larger, and the phonon thermal 
conductivity can be killed off by cutting away the 
long-wavelength phonon paths along the c-axis. 



D. More realistic estimates 

To further validate our results, we can give an 
estimate of the overall zT (up to now, we have 
studied the electronic part only). Considering the 
Wiedemann-Franz relation, the figure of merit can 
be rewritten as: zT= ctTS 2 /k= S 2 /Lq being Lq 
the Lorenz number with a value for free electrons 
L = 2.45 x 10~ 8 (V/K) 2 . At high enough tempera- 
tures, the lattice thermal conductivity term typically 
decreases as 1/T (with the electronic term roughly 
temperature independent ). This type of behavior 
can be observed in Ref. 23. Taking the values at 
the transition temperature for the thermal conduc- 
tivity from Ref. [23|, the electrical conductivity from 



Ref. [H, and using our calculated values for the ther- 
mopower, we can obtain a more realistic estimation 
of the zT values for each functional at each doping 
level. The behavior of the estimated zT with tem- 
perature for the various doping levels and functionals 
used can be seen in Fig. [6] It can be observed that 
the shape of the curves for each functional remains 
roughly the same than in Fig. [5] The calculated 
electronic-only factor ctT/k is not far from constant 
when the real band structure of the material is intro- 
duced in the calculations. The overall values of zT 
are reduced with respect to those obtained only with 
the electronic contribution in approximately one or- 
der of magnitude due to the phonon thermal con- 
ductivity. zT ~ 0.1 is predicted from our calcula- 
tions at high temperature if the proper oxygen con- 
tent is chosen. In any case, the growth in the form 
of thin films with the appropriate oxygen content 
would increase the electrical conductivity and reduce 
the phonon thermal conductivity. Our results show 
that, in that case, the system could become promis- 
ing for further research in terms of its thermoelectric 
properties, trying to find other oxides with perfor- 
mance comparable to that of layered cobalt ates. 



IV. SUMMARY 

Our ab initio calculations for the compound indi- 
cate some promising features of hole-doped La2Ni04 
as a possible oxide thermoelectric compound. If the 
thermal conductivity can be reduced by nanostruc- 
turing, e.g. in the form of thin films, the system 
could show enhanced thermoelectric performance at 
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low hole-doping levels, attainable by the appropri- 
ate control of the apical oxygen content. Our cal- 
culations show that a region with relatively large 
Seebeck coefficient exists in this compound at small 
doping levels, within the realistic AF description. A 
careful experimental study needs to be performed in 
this respect controlling the oxygen content, and also 
making thin films with the appropriate oxygen com- 
position. Hole-doping will increase the conductivity, 
as the thin film geometry does, which together with 
the reduction in thermal conductivity could leave 
room for an improvement of the thermoelectric fig- 
ure of merit of this layered nickelate. In terms of 
the electronic structure, band engineering could be 
used to explore the possibility of using the large ther- 



mopower expected from the d z i band being closer to 
the Fermi level. 
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